Field enhancement factor and field emission from 
a hemi-ellipsoidal metallic needle 
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We present an exact solution for the electrostatic field between a metallic hemi-ellipsoidal needle 
on a plate (as a cathode) and a flat anode. The basic idea is to replace the cathode by a linearly 
charged thread in a uniform electric field and to use a set of "image" charges to reproduce the anode. 
We calculate the field enhancement factor on the needle surface and ponderomotive force acting on 
the needle. Using the Fowler-Nordheim theory we obtain an exact analytical formula for the total 
current. 
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I. INTRODUCTION 

Application of various one-dimensional nanostructure 
materials as field emission sources has attracted exten- 
sive scientific efforts. Elongated structures are suit- 
able for achieving high ficld-emission-current density at 
a low electric field because of their high aspect ratio. 
Area of its application includes a wide range of field- 
emission-based devices such as flat-panel displays, elec- 
tron microscopes, vacuum microwave amplifiers. X-ray 
tube sources, cathode-ray lamps, nanolithography sys- 
tems, etc. 

Since the discovery of carbon nanotubes (CNT) [1-3] 
and experimental observations of their remarkable field 
emission characteristics [4-6] , significant efforts have been 
devoted to the application of using CNT for electron 
sources. Recently, field emission from metals [7], metal 
oxides [8-11], metal carbides [12], and other elongated 
nanostructures have also been explored. It is now pos- 
sible to control the diameter, height, radius of curvature 
of the tip, and basic form of emitters during growth. 
Elongated structures of different shapes such as nan- 
otubes, nanocones, nanofibers, nanowires, nanoneedles, 
and nanorods have been successfully grown [13-17]. 

The current emitted by elongated nanostructures de- 
pends on the distribution of field strength at its top. A 
great deal of attention of researchers has been devoted to 
the apex field enhancement factor. Various approxima- 
tions were used to evaluate the field amplification. Typi- 
cally, the field enhancement factor is calculated from the 
slope of the FowlerNordheim plot. Under such considera- 
tion one can figure out only the field enhancement factor 
averaged over the entire emitting surface. Various theo- 
retical approaches based on both analytical solutions and 
numerical simulations have been reported. [18-24] Analyt- 
ical solution was obtained for a hemisphere on a plane, 
a floating sphere at the emitter-plane; a hemi-ellipsoid 
on a plane [18], and a sphere-on-cone model [19]. Usu- 



ally only the apex field enhancement factor is calculated 
when dealing with a hemi-ellipsoid on a plane [18]. In- 
fluence of the anode-cathode distance within a floating 
sphere model was investigated in [20]. The model of a 
hemisphere on a post for CNT emitters was used in nu- 
merically simulations and reported in many papers [21- 
24]. Calculation difficulties in these numerical methods 
arise due to the large nanotube aspect ratio and very long 
distance between cathode and anode in comparison with 
emitter height. Usually, these numerical results were gen- 
eralized and simple fitting formulas of field enhancement 
factor for individual nanotube [21, 22, 25, 26], for nan- 
otube in space between parallel cathode and anode planes 
[27-30], and for a nanotube surrounded by neighboring 
nanotubes with a screening effect [31-35] were suggested. 

The main problem for such algebraic fitting formulas 
is the lack of a definitive proof of their accuracy. Only 
a few papers considered forces acting on nanoemitters 
under electric field [32, 36-38]. Thus far, there is no ana- 
lytical formula which provides a good approximation to 
the total current generated by the nanoscale field emit- 
ter. In this paper, we consider the electric field strength, 
field enhancement factor, ponderomotive forces, and to- 
tal current of a metallic elliptical needle in the form of 
hemi-ellipsoid in the presence of a flat anode. We pro- 
pose an exact analytical method, exploiting a very simple 
idea. The basic idea is to replace the cathode by a linear 
charge distribution in a uniform electric field and to use 
a set of "image" charges to reproduce the anode. 
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II. BASIC IDEA OF ANALYTICAL METHOD 



Consider a prolate metallic spheroid in a uniform elec- 
tric field. We can replace the spheroid by a linearly 
charged thread. The thread is a green line on Figure 
1 and the linear charge distribution is represented by a 
red line. The length of a thread is 2h. The electrostatic 
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FIG. 1: Linearly charged thread in a uniform electric field 
along z. 



potential produced by the charged thread is: 

h 

ip{r, z) =- 



1 



'dz' 



(1) 



where Sq is the electric permittivity, (r, z) denotes the 
in-plane radial and z coordinates, tz is the linear charge 
density at point (0,z), h is half of the thread length. 
The solution is independent of the azimuthal angle. We 
assume the coefficient of linear charge density, r to be 
positive. 

After integration we have 
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The total potential is 

$(r, z) = ip{r, z) + EqZ 
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where Eq is the applied macroscopic field or far field. The 
shape of the metallic spheroid is given by the solution to 
the equation 



(y9(r, z) + Eqz — 0. 



(4) 



Using new variables — \J (z — h)'^ + and rh = 

{z + hY + and the equality z = ° , we can 
rewrite Eq. (4) in the form 
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where C — 47r£o-^- From this equation we obtain Va + 
rfc = const, which means the cross-section of the spheroid 
is truly an ellipse. Points (0, —h) and (0, h) are the focii 
of the ellipse. Therefore and rf, are distances between 
(r, z) and the two focii. 



III. GEOMETRY OF THE ELLIPSOIDAL TIP 
OF THE METALLIC NEEDLE 



Using a dimensionless parameter, the eccentricity 

Ta + n 

we transform Eq. (5) into 



C = In 
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(6) 



(7) 



If ^ is close to 1, the ellipse becomes elongated. If ^ ^ 
the ellipse turns into a circle. Therefore, by changing the 
coefficient of linear charge density, t we may modify the 
shape of the ellipse. 

We can also adjust other geometrical parameters of 
the ellipse: the length of semi-major axis or height H\ 
the length of semi-minor axis or base radius, i? at z = 0; 
and radius of curvature, p at point (0, H) (see Fig. 1). 
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IV. ELECTRIC FIELD AND FIELD 
ENHANCEMENT FACTOR ON THE SURFACE 
OF THE METALLIC SPHEROID 

Taking the gradient of Eq. (3) and using Eq. (4) we can 
calculate components of the electric field on the surface 
of the metallic spheroid: 
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Thus the modulus of the electric field is 
Eo 4h^{n-ra) 
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(11) 



Eq. (9-11) allow determining the electric field strength 
on the surface of the half ellipsoid at an arbitrary point. 
The field enhancement factor at the apex of the ellipsoid 
follows from Eq. (11) 
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Analytical expressions for field strength on the z axis 
and for field enhancement factor on the tip of the half 
ellipsoid obtained previously [18, 39-41] are in agreement 
with our result. Here, by taking gradient of Eq. (3) we 
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FIG. 2: Comparison of the exact formula (12) with approx- 
imate formula (13) for the field enhancement factor versus 
H/p. 



can obtain the field strength at any point we desire. In 
the limit ^ ^ we have a metallic half sphere and the 
field enhancement factor = 3. If ^ ^ 1 then for the 
elongated metallic needle, we have 
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In Fig. 2 we see that the approximate Eq. (13) has a 
vertical asymptotic line at H/p = e^/4 « 1.85. The 
exact Eq. (12) gives 3 for the half sphere when H/p 1. 
However, when H/p > 100, the relative error becomes 
less than 1.7%. 




FIG. 3: Geometry for the calculation of ponderomotive force 
acting on the belt between ra = A and ra = B. 



Using ra — H — zS^ wc rewrite 
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where z varies from za to Zb (see Fig. 3). It gives us 
the value of net detaching force acting on the surface of 
the spheroid integrated over the surface between planes 
z — ZA and z — zb- The total detaching force acting on 
the ellipsoidal needle is 



PONDEROMOTIVE FORCES 
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Let us estimate the electric force acting on the surface 
of ellipsoid. The electrostatic force acting on the elemen- 
tary area, s of the external surface is given by 



F 



E fids. 



(14) 



where n is a vector normal to the surface. 

We can calculate the force acting on the spheroid sur- 
face between circles Va = A and Ta = B (see Fig. 3). It 
is obvious that r-component of force on that surface is 
equal to zero. For the z-component, after routine opera- 
tions we obtain 



dF, = ^-^E^n.ds = ^E^dr^ = ^E^—dva. (15) 
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We can integrate the above equation analytically for 
changing from AtoB. We have 

FM\B) = ^-^^^^^^^ 
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In Fig. 4a we show the relative detaching force 
Fz{H, z) / Ftotal as a function of coordinate z, where 
Fz(H, z') is the force acting on the surface between plane 
z = z' and the tip of the spheroid (z = H). Distribu- 
tion depends only on the parameter H/p. The major 
part of the detaching force is concentrated near the tip 
when H/p is large. Fig. 4b shows the total force isolines 
on the {p,H/p) plane. When we chose logarithmic scale 
for p and H/p with logarithmic steps we obtain a set of 
nearly straight isolines with equal distances in the plot. 



VI. FIELD EMISSION FROM INDIVIDUAL 
NEEDLE 

The field emission from a metallic surface is described 
by the Fowler-Nordheim theory. According to this theory 



exp — 



(19) 



where j denotes the emission current density in A • cm^^, 
E is local electric field at the emitting surface in V • 
cm~^, (b is work function in eV, and the first and second 
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FIG. 4: Distribution of force and total force isolines. (a) Distribution of relative force over the axis of needle, (b) Isolines of 
total force on the plane. 




p, nm 

FIG. 5: Distribution of emission current and total emission current isolines for work function (j) — 4.8 eV (a) Distribution of 
relative emission current over O. (b) Total emission current isolines on the (/3i5o,p) plane. 



FowlerNordheim constants are Ci = 1.56 x 10 ^[A • eV • 
y-^], C2 = 6.83 X 10^[V • eV"^/^ • cm"!] respectively. 

To integrate Eq. (21) over the surface, we express 
the elementary area as the following ds = 2'Kydl = 
2TTr a sm Or a dO — —2Trr^dX, where A = cos 6*, 9 is an- 
gle between axis z and as shown in Fig. 1. Ellipse 
equation in polar coordinates is 

I 



where A G (—^,1) gives half of the ellipse which we are 
interested. In the case of half sphere we have 



(3 = 3, 3EoX, p = H. (21) 



The current emitted from the surface bounded by G 
(0,6) is 



j(e) = 2V 



A^ exp 
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where C-< 
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and El (a;) = / 



r cxp( — 



dt is the ex- ponential integral. Due to small /3 for the half sphere we 
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need to use very strong electrical field to produce slightly 
visible current in experiment. In the general case, the ex- 
act formula for emission current should be written as 
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where for Tq we use relation (20). Unfortunately it is not 
possible to integrate (23) in the general case, but we can 
obtain approximate formula for the elongated ellipsoidal 
needle. In this case the emission takes place mostly from 
the tip of the needle where the field is max;imal due to 
the curvature being the largest there. 
Due to p <^ H we have 

P 



1 + A' 



(24) 



For the evaluation of (23) we assume H — Va ~ H , 2iJ — 
Ta « 2H. Then the integral can be expressed in terms of 
the El-function and we get 



j(e) = 



2Trp^C,il3Eof 



El 



(25) 



By comparing with accurate numerical integration of (23) 
we find that formula (25) is accurate up to the first four 
digits for H/p > 100. So formula (25) is a good approxi- 
mation for investigating the current emitted from an area 
of surface depending on angle 0. The total current with 
similar accuracy can be written approximately as 
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The emission current has much sharper distribution on 
the tip than the detaching force due to exponential de- 
pendence of Fowler-Nordheim formula (19). In Fig. 5a 
we plot the relative emission current, J {Q) / J total as a 
function of (angle between Va and axis z), where J(0) 
is emission from the tip on the surface of the ellipsoidal 
needle. From (26) we see that total emission current de- 
pends only on the area of the hemisphere, 27r/9^, work 
function, ((>, and the electric field on the tip, (iEo. So 
in Fig. 5b we show the total emission current isolincs 
within a reasonable range (from IQ-* to 10"* nA) with 
logarithmic steps on the {/3Eo,p) plane, where the work 
fimction (f> = 4.8 eV is assumed. If the electric field is 
Eq « 1..5 V//xm, then we have to use an enhancement 
factor (3 = 10^. .10* to get an appreciable current. 



fluence on the emission current and the detaching force. 
The basic idea of calculation is to replace the cathode by 
a linearly charged thread in a uniform electric field and 
to use a set of "image" charges to reproduce the anode 
as shown in Fig. 6. We put infinite set of "images" of 
the linearly charged thread with the same spacing, 2L 
{L > H). It is clear that planes z = Q and z = L are 
planes of symmetry for distributed charges and will have 
potentials on cathode and V = EqL on anode. With 
the same r we will get different geometry of cathode due 
to additional potential of "images" . We assume that the 
new form of the elongated needle will be approximately 
described by ellipsoid especially near the tip. On the sur- 
face of the thin ellipsoid, we have r <^h and z < H < L. 
So, we can describe the image potential ipi as a function 
of z. 
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where C* = iTreoEo/r*, L = r]H, rj > I. The thread 
potential is 
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Near the tip z ^ H wc have the following equation for 
describing the shape of the ellipsoid: 

ip{r, z)+Eoz+ipi{z) « ip{r, z)+EQZ+^^pi{H) = 0. (29) 



Then in analogy to (5) we write 
4/i , frn.+rb + 2h 
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where 
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So we have 



C* =ln( l±i ) _2^-P. 



Analogous to (9) and (10) we have 



(31) 



(32) 



VII. SIMULATION OF PLANAR ANODE BY 

"IMAGE" CHARGES E;, 

The presence of a flat anode placed at a distances com- 
parable with the length of the nano needle has strong in- 
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FIG. 6: Infinite set of "image" charges for the simulation of a planar anode. 
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FIG. 7: Isolines for the field enhancement factor f3 on the 
{L/H,H/p) plane. 



where 



C* rany/ira + rf,)^ -4/l2 



(34) 
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If |(P - Q)/(/3 • C*)| < 0.01 then we may neglect {P 
Q)/C* w in (33). We get 



PiL/H,H/p) 
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Finally for the electrical field with the assumption 
[P - Q)/C* « formulas (33) and (34) are completely 
analogous to (9) and (10). Therefore for the force and 
emission current we may use the above derived formulas. 
Instead of (12) we should use (31), (32), and (36). So, the 
parameter rj = L / H in the force and emission current for- 
mulas (17), (18), (25), and (26) should be included only 
through the field enhancement factor /3 — (3{L/H,H/p). 

In Fig. 7 we present isolines for the field enhancement 
factor P on the {L/H,H/p) plane. The dashed line di- 
vides the graph into the upper-right and lower-left side. 
On the upper-right side our approximation has less than 
1% error, |(P - Q)/(/3 • C*)| < 0.01. 

What distance between the anode and cathode is large 
enough to assert that the elongated metallic spheroid is 
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FIG. 8: Emission current as a function of the anode-cathode 
distance parameter 77 = L/H for constant electric field Eq = 
V/L (red line) or Eq = V/{L - H) (blue line). Radius of 
curvature is fixed at p = 5 nm. For solid curves, H/p — 5800 
and for dash curves, H/p — 10000. For red and blue solid 
curves Eq = 5V//im, for red and blue dash curves Eo = 
3V//im. 



placed in a uniform electric field? In experiment we may 
measure distance L between the anode and cathode and 
the distance, L — H between the anode and the needle 
apex. We can determine the applied electric field both by 
Eq — V/L and by E^ — V/{L — H). In our calculations 
we move the anode and increase the voltage V so as to 
keep Eq or Eq immutable. It is clear that for a large 
distance (in the limit L — s- 00) the difference between 
definitions of Eq and E[^ disappears. Fig. 8 shows the 
total current emitted from the needle versus the anode- 
cathode distance. 



We set for blue lines Eq = const and for red lines Eq = 
const. With Eq = Eq — const the total currents have 
the same limit for 7/ = L/H 00. But currents with 
Eq = const tend to approach the limit much slower than 
the corresponding currents with Eq — const. If we define 
the applied electric field as Eq and the anode-cathode 
distance is ten times more than the needle height, then 
we can neglect the influence of the anode and assume 
that the metallic needle is placed in a uniform field. In 
contrast, if the applied electric field is defined as E'q^ then 
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the anode-cathode distance must exceed one thousand 
times of the needle height for the above statement to be 
vaHd. 



VIII. CONCLUSION 

In this work we theoretically investigate the field emis- 
sion from an elongated prolate spheroid metallic needle 
in diode configuration between a fiat anode and cathode. 



Exact analytical formulas of the electrical field, field en- 
hancement factor, ponderomotive force, and field emis- 
sion current are found. Applied voltage, height of the 
needle, radius of curvature on its top, and the work func- 
tion are the parameters at our disposal. The field en- 
hancement factor, total force and emission c;urrcnt, as 
well as their distributions on the top of the needle for 
a wide range of parameters, have been calculated and 
analyzed. 
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